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I.  INTRODUCTION 


Using  simple  friction  models  in  digital  computer  simulations  sometimes  leads  to  less  than 
satisfactory  results.  Consequently,  two  friction  models  were  developed,  one  being  a 
modification  of  the  other.  These  models  essentially  combine  two  simple  friction  models. 

The  two  models  are  compared  with  a third  which  is  often  used  for  friction  modeling.  For 
testing,  the  models  were  embedded  in  a second  order  spring-mass-damper  system. 
Comparisons  are  made  between  the  three  systems  for  different  system  gains,  friction  values 
and  damping  using  sine  wave  and  step  inputs.  Only  the  RK2  integration  algorithm  was  used 
because  of  its  predominant  usage  in  the  author's  simulations.  The  test  program  was  written  in 
the  Advanced  Continuous  Simulation  Language  (ACSL).1 

2.  MATH  MODELS 

The  friction  models  are  incorporated  into  a second  order  spring-mass-damper  system. 
Similar  system  models  are  often  used  to  represent  missile  seeker  platform  dynamics.  Figure  1 is 
a block  diagram  of  the  test  system.  The  closed  loop  equation  without  the  friction  nonlinearity 
is 

Jp_  _ K 

TM3  S2  + B*  K • S + K • CMPL 

where  t h is  the  platform  angle  ( radians);  TMj  is  the  motor  torque  into  the  system  (ft-lb);  B is  the 
viscous  damping  (ft-lb/ rad  /sec);  CMPL  is  the  spring  constant  (ft-lb/ rad);  and  K is  the  inverse 
of  the  platform  inertia  (lb-sec’-ft)1.  The  values  for  TMt  and  CMPL  were  chosen  to  keep  t/r  well 
below  tt/2  rad.  Initially,  TM,  is  5 ft-lb  and  CM  PL  is  12  ft-lb/ rad.  Bis  made  0.1  ft-lb/ rad,  sec, 
which  gives  a closed  loop  damping  ratio  of  about  0.08  for  a K of  30  (lb-sec’-ft)'1 . For  all  three 
friction  models  a value  of  2 ft-lb  is  initially  used. 


I f.  I Mitchell  and  Joseph  S.  Gaulhier.  Advanced  Continuous  Simulation  Language!  ACSL)  User  Guide  Reference  Manual. 
Mitchell  and  Gauthier  Associates.  1975. 
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Ollen  tiiclion  modeled  as  shown  in  Figure  ' In  M.S1  code  the  lnetion  model  is 
Iik  *S(iNi-(i'i  I his  model  w ill  he  i elei  i ed  to  .is  model  \ \ diaw  hack  ol  this  model  is  that  it 
ptoduees  .iitdieial  torques  when  there  is  no  jnmh.il  motion  Model  A is  used  loi  test 
eompaiison  purposes  I he  two  lnetion  models  designed  to  improve  upon  this  model  are 
shown  m h\;urt‘\  > and  4 I lies  will  he  ealled  models  It  and  (.'.  respeeti\el\ 


In  M’Sl  eode.  lnetion  model  If  is  meoi  pointed  into  the  expression  lot  net  torque  as 
lollows 


tnet  - RSW 


vRMIN,DEAD(_TFR1'  TFR1'TM3-B’‘^  ~ CMPL*^ 


TNET  + SIGN(TFR2'  - **>> 

W hen  li.^the  gimbal  angulat  rate,  is  less  than  the  value  ol  Rmis,  the  gimbal  friction  ean  he 
modeled  as  a dead  /one  since  no  {timbal  motion  will  occur  until  I m i exceeds  a value  equivalent 
to  the  trietional  resistance  torque.  1 i ki  Once  the  jumbal  is  in  motion  the  {timbal  friction  will 
produce  an  effective  torque.  Iik-,  that  always  opposes  jumbal  motion.  However,  false  motion 
can  occur  due  to  at  least  two  things:  (I)  i//will  not  necessarily  be  smaller  than  Rmis  (the 
incremental  change  in  lA^ovet  one  integration  step  can  be  larger  than  Rmis)  so  that  the  dead 
/one  model  may  not  switch  back  in.  (I)  If  the  dead  /one  model  does  not  switch  back  in.  the 
trietional  torque  ean  not  only  oppose  motion,  hut  can  also  cause  a reversal  ol  gimbal  motion, 
w Inch  is  not  physically  correct  So  an  appropriate  value  must  be  chosen  for  Rmis  for  model  H 
to  operate  properly  This  is  an  undesirable  feature  that  was  eliminated  in  the  design  of 
model 

As  shown  in  Figure  4.  model  is  comprised  ol  model  B plus  logic  that  zeroes  the  gimbal 
angular  rate,  i A.  whenever  both  the  net  torque  minus  friction  is  less  than  the  friction  and  tA 
changes  sign  It  should  he  noted  that  the  dead  /one  model  uses  1 1 m.  and  the  sliding  friction 
model  uses  1 1 « • I his  allow  s foi  friction  models  that  have  different  v allies  loi  static  and  sliding 
friction  I he  logic  for  switching  iA  to  zero,  as  implemented  in  / ignrr  4.  forms  an  algebraic 
loop  I hus.  model  C was  implemented  as  an  Af'Sl  macro  in  which  the  loop  is  broken  at  the 
desired  point  in  the  system  model 


I he  macro  code  and  a How  diagram  ol  the  macro  are  shown  m Figures  5 and  b.  respectivelv 
In  the  macro.  iA  (macro  local  variable  VO)  is  not  actually  set  to  zero  except  w lien  I OTORQ  is 
I Rif  ias  discussed  below ) <A  is  calculated  hv  adding  <An(  YHO  H and  iAi>0  OO  I P)  as  show  n 
in  Figure  y Setting  i h to  zero  through  biasing  ifr,,  vv  it h iin\ ( is  set  to  -i/m)  keeps  tA  from  hav  mg 
a tump  in  its  value  whenever  1 OIORQ  goes  Irom  I Rl  I to  I Al  SI  . A pimp  in  ii<  would 
otherwise  occur  since  i//|>.  being  a state,  cannot  he  set  hack  to  zero  l O I ORQ  becomes  1 Rl  l 
whenever  Ki  i.  the  net  torque  minus  friction,  i'  less  than  Iik:  and  iA  has  switched  sign. 


S 


When  the  A1  SI  I xecutivc  in  in  the  INI  I IAI  sectuin(/./IC  M isTRUF)theexpressionfor 
I O I ORQ  is  bypa  sed  since  Ism  will  not vet  be  defined  When  Z7.ICF1.  is  FALSE,  LO  TORQ 
is  calculated  as  discussed  above  II  I O I ORQ  is  1 R 1' I . i/i  is  set  to  zero.  II  this  integration  step 
is  the  first  minor  step  ol  a multi-step  integration  method  (//.I  S I = 1.0),  then  i j/[,  is  set  to  -<Ai>. 
<l/\t  is  switched  only  once  lor  each  multi-step  integration 

When  1.01  ORQ  is  FAFSF  or  Z/.IC'H  is  1RUF.ZZFS1  is  checked  for  a value  of  1 .0.  If 
ZZFS1  is  1.0,  then  YDl  , the  old  value  ol  t h.  is  set  to  i J/. 

Ts.ii.  the  net  torque  including  friction,  and  t[n>  are  calculated  outside  of  the 
I’ ROC  F DU  R Al.  block  the  rest  ol  the  code  is  in.  The  expressions  for  Tm?  and  4i\>  are  not  in  the 
block  in  order  to  allow  them  to  be  sorted  by  the  ACSI.  compiler  along  with  the  expression  for 
Fsi  i.  In  the  FOR  IRAN  compilation  of  the  test  program  using  the  FRICTN  macro,  the  code 
in  the  PROCEDURAL  block  came  first  followed  by  the  expressions  for  Tm  i,  Tm  i and 
respectively. 

«i»/>  and  YDL  are  only  reset  once  per  multi-step  integration  to  keep  Tsh  from  having 
erroneous  values.  An  example  of  how  problems  occur  when  ij/i  and  YDl.  are  reset  more  than 
once  is  shown  below.  For  RK2  as  implemented  in  ACSL’s  FORTRAN  library. 

( t + h)  = ^(t)  + 

+ 3*'PD(t  + J h>VD  (t)  f j h (t)) 


where 

h is  the  major  integration  step  si/e 
$r>  = K*  T\ei 

Define  i /rm  to  be 


*D1  = (t  + Th'  (t)  + Ih  (t) 


Suppose  LOTORQ  goes  TRUE  on  the  intermediate  integration  step.  Then  i/rm  will  be  zero 
since  Tm  i and  i h are  zero.  Suppose  tim  is  set  to  -tim  and  YDL.  is  set  to  i i>  on  this  intermediate 
step  On  the  final  evaluation  for  ij/i>  (t  + h)  the  RK2  algorithm  reduces  to 
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<I>D  (t  + h)  = i^D  (t) 


if>D  (t,  i|^n  (t)) 


D 


For  LOTORQ  to  have  gone  1 RUF.  on  the  intermediate  step, 
*D  (t)  + Ih  (t)> 


the  expression  for  t/in  (t  + h ).  had  to  be  of  opposite  sign  from  i/»n(t),  since 


* = + *D 

(Assume  ij/!>  = 0.0  before  I.OTORQ  went  TRUE). 

Assuming  iji » (t  + t h).  is  negative  and  t)  is  positive  then  the  expression 

(t,  $»_  (t)) 


>PD  (t  + h)  = ipD  (t)  + 
can  conceivably  be  positive  or  negative. 

Suppose  it  is  negative.  Then,  on  the  next  derivative  calculation  for  ij/D  (t  + h). 


tf)  (t  + h)  = (t  + h)  + tj»' 

- *D  (t  + h)  - <j>D  (t  + |h  ) 

? 0.0 

Since  YDL  is  zero  when  I.OTORQ  is  calculated,  I.OTORQ  will  be  FALSE.  ThusTsM  and  i^(t 
+ h)  will  have  nonzero  values  and  i^o  (t-f2h)^  (t-fh).  So  ij>  is  not  zero  as  it  should  be.  and  &> 

is  continuing  to  be  integrated. 


To  make  the  system  behave  properly.  i^Sand  YDL  must  not  be  reset  on  ari  intermediate  step 

since  the  value  of  i^>  (t  + h)  does  not  depend  only  on  im.  It  is  to  be  noted  that  either  the 
• • • 
resetting  of  t/n>  or  YDL  would  have  caused  problems.  If  i/n  and  t/n>  had  not  been  reset 

(t  + h)  = (t  + h)  + iji' 

= (t  + h)  + 0.0 

^ 0.  usually 
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Since  > I >1  is  i//(()  I 0i,m  whemintt)  O.Otwhen  I OIORQ  is  calculated,  then  IOIORQ  will 
he  1 Rl'l  if 


0D  (t)  * 0D  (t  + h)  < 0 

and0tt  > h)  will  be  reset  liont  0i>  tt  + hftozero  1 m i ami  0n(t  1 h)  will  consequently  he  zero, 
ami  0D  (t  + h)  will  remain  constant.  II 

•I’d  (t)  * 0D  (t  + h)  >0 

then 

& (t  + h)  ft  0, 

ami  the  system  will  continue  to  integrate  0d  .is  it  should 

3.  TEST  PROGRAM  AND  RESULTS 

Figure  7 is  a listing  of  the  test  program  In  the  1>I  RIVA I I VI  section  seeker  platform  angles 
SI  (0),  SIP  and  SIPP  are  calculated  lor  friction  models  If  and  A.  respectively.  1M1  is  the 
input  motor  torque  for  all  three  models  If  SINS  IP  is  I Rill  I isa  sine  wave;  otherwise.  1 M< 
is  a step  function.  I m i and  1m  r arc  the  net  torque  minus  friction  for  models  C and  II, 
respectively.  I m i and  I m h are  the  net  torques  with  friction  for  models  Rand  A,  respectively. 

Parameters  and  initial  conditions  are  defined  in  the  INI  1 1 A I section  through  I'ONS  IAN  I 
statements.  I he  TIN  1 1 RVAI  statement  specifies  the  data  recording  interval.  1 lie  RK2 
integration  algorithm  is  specified  by  setting  I A I li  to  4 in  the  Al  liORI  I MM  statement.  I he 
M AX  1 1 RVAI  statement  overrides  the  NS  I l l’S  statement  to  specify  a 2.5  msec  integration 
step  size.  WP  is  the  frequenev  of  the  l M,  sine  wave  in  rad  sec  and  W is  WP  in  11/ 

ATS  I run  time  commands  are  used  for  altering  parameters  and  initial  conditions  on  states; 
running  the  simulation;  and  specifying  the  output.  I he  simulation  was  run  intcraetivclv  on  a 
TV  RPR  74  using  run  time  commands  entered  on  a lektroni.x  console.  Hard  copv  plots  were 
generated  along  with  line  printer  listings  Run  time  commands  that  were  used  evert  session 
were  put  on  a temporary  mass  storage  file  that  could  be  attached  as  a local  file  to  the  user's 
interactive  terminal 

Figure  <s’  is  a list  of  interactive  session  setup  and  ACS!  run  time  commands.  After  logging 
onto  the  l V RFR  74  IN  1 ERl'OM  system  via  the  1 ektronix  terminal,  the  following  setup 
statements  are  entered: 

1 1 


r 


CONN  I C l.  Ol  I I’ll  I 

specifics  the  system  output  tile  dotau It  name  OH  I I’ll  I.  will  he  displayed  on  the  lektronix 
screen 

Eli  . 170 

extends  the  execution  time  per  statement  entered  to  170  octal  sec. 

ATTACH.  INPUT.  TF.K.  ID  = DDXXXII 

attaches  the  highest  cycle  (2)  ol  I I K to  the  I eklronix.  IN  I’ll  I is  the  default  file  name  for  input 
to  local  programs  executing  in  the  system. 

ATTACH.  I GOB.  TEK.  ID  = DDXXXH.  CY  = I 

attaches  the  simulation  absolute  binary  file  which  was  previously  compiled  in  a batch  job  via 
the  system  control  cards  listed  in  Figure 

After  this  initial  preparation.  EGOB  is  executed  by  the  statement 

l GOB 

and  the  simulation  reads  the  ACSl.  run  time  commands  on  file  INPUT,  the  contents  of  which 
are  shown  in  Figure  Vb: 

SET  PRN  =9 

tells  the  simulation  that  the  line  printer  data  will  he  placed  on  the  PRINT  file  (logical  unit 
number  9)  instead  of  the  default  file  name  OU  I PU  T.  After  termination  of  the  simulation  the 
PRINT  file  is  batched  to  a line  printer. 

SET  PRN  PET  = .1  ..  CAI  I’l  l = .T„  1 1 1 CPE  = .T. 

replaces  the  default  printer  plot  routines  with  the  special  lektronix  plot  routines  ami  ! 1 1 CPI 
set  TRUE  causes  titles  to  be  placed  on  the  lektronix  plots. 

SEE  HU  E = ’ERICTN  IIS  1ER-* 
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pecilics  llu'  lust  two  woids  i>l  llu'  plot  title 


SI  I m MP  I 

causes  the  Iasi  \ allies  ol  all  the  \ai  tables  to  he  wi  men  to  the  lile  with  logical  unit  nuinhei  PR  N 

SAN  I H ASI  1 

‘saxes'  the  initial  parameters  amt  I V ‘s  lot  latei  use 

I’RtK'l  D (j() 

SI  AR  I 
PRIM  I ‘All* 

INI) 

ileli  ties  a I’RtK'l  I > that  stmts  execution  ol  the  simulation  ami  causes  the  \ at  tables  specified  in 
the  I’KI  PAR  command  to  he  recorded  even  ( IN  I seconds  onto  the  lile  with  logical  unit 
number  I’RN 


PRI  PAR  I . SI.  SID.  SIP.  SU’D.  SIPP,  SII’PD.  I M l. 

I Ml.  INI  12.  INI  1 I 

specifics  the  variables  whose  values  will  be  output  on  a line  printer. 

SI  I t Ml)  DIS 

indicates  to  the  simulation  that  i un  time  commands  will  now  be  input  from  the  DIS  PI  AY  file 
whose  logical  unit  number  is  DIS  DIS  is  h\  default  h which  is  also  the  numhet  lot  the 
OH  I PlH  lile  I luis,  i un  tune  commands  are  now  expected  to  come  from  the  lektronix 
terminal 


SI  I I IK  I I4.IIR2  14.0  20.  K It).  MAXI  0.01 5 

Although  MAXI  was  set  to  0 015  the  integration  step  si/e  was  0.010  see  due  to  (INI  being 
0 010  Instead  ol  starting  w ith  a typical  spring -mass -damper  system,  the  pat  ameters  have  been 
changed  to  reflect  how  they  would  be  relative  to  each  othci  lot  a missile  seeket  plattorm 
Namely,  the  frictional  toriptes  arc  highet  than  ton|iies  due  to  the  spt  mg  constant  and  v iscotis 
damping;  K the  reciprocal  ol  the  platform  inertia,  is  the  system  gain  and  is  at  bit  rat  tlx  set  to  10, 

I l 


;nul  (i  in  made  20.  It-lbs  in  ordci  lor  the  platlorm  lorquci  to  overcome  the  friction.  MAX  1 was 
changed  to  reflect  the  lower  value  ot  the  undamped  natural  frequency. 

SI  1 111  I 1 1 3 1 ‘8  .HIM  7<T 

fills  the  third  word  ol  the  plot  title  with  the  dale. 

C.O 

starts  the  simulation  through  PROC'FI)  (i(). 

plot  ‘xnr  = 2..  si.  sin 

limits  the  plot  abscissa  to  2 sec  and  causes  the  plotting  of  SI  and  SID  (t/r  and  t/d  for  the  system 
with  model  0.  Figure  10  is  the  resulting  I cktronix  hard  copy.  SI  limits  to  a steady  state  value 
and  SI  I > stays  at  zero.  SIP  and  SU’D  (i //  and  < //)  for  the  system  with  model  B and  SIPP  and 
S1PPD  ti h and  >ii)  for  model  A arc  plotted  in  Figures  1 1 and  12.  respectively.  Notice  SIP  and 
SIPP  do  not  reach  a steady  state  value  and  that  SU’D  and  SII’I’D  are  chattering  when  they 
physically  should  not  be.  Figure  1.1  compares  SI.  SIP  and  SIPP.  SIPP  leads  SIP  and  SI 
primarily  due  to  the  initial  spurious  oscillations  in  SIPPD  before  IMi  rose  above  zero.  Since 
RMN  for  model  Bis  set  to  It)  in  the  INITIAL  section  ofthe  program,  spurious  oscillations 
arc  to  be  expected.  RMN  was  reset  to  0.0 1 5 rad  sec.  and  the  resulting  SIP  and  SU’D  arc  shown 
in  Figure  !4a.  Only  an  improper  step  in  SU’D  shows  up.  It  is  due  to  the  jump  in  TN  FTP  caused 
by  SU’D  changing  sign  and  still  being  greater  than  RMN  in  absolute  value.  Eventually  SU’D 
gets  caught  in  the  ±RMN  interval  and  model  B reverts  to  a dead  zone  model.  Figure  14b 
compares  the  angles  of  models  B and  C for  RMN  = 0.015. 

Resetting  RMN  to  10  '"and  setting  SINS  IP  to  TRUE  (the  torquer  motor  input  is  a 1.0  Hz 
sine  wave)  results  in  the  plots  shown  in  Figures  1.1- 18.  Note  that  although  SU’D  and  SII’I’D  do 
not  stay  zero  when  they  should,  the  comparison  of  SI.  SIP  and  SIPPD  in  Figure  IS  is  not  as 
bad  as  for  the  step  inputs.  The  systems  with  models  A and  B showed  improved  sine  wave 
response  with  a decrease  in  step  size,  but  their  step  responses  still  did  not  steady  out  when  the 
step  size  was  lowered. 

RMN  was  set  to  0.5  tad,  sec  and  the  resulting  SIP  and  SU’D  for  model  B arc  shown  in 
Figure  IV.  Figure  20  compares  SI  and  SIP. 

With  a call  to  "RFSTOR  ‘BASF”'  the  original  parameter  and  l.t\  specifications  in  the 
INI  MAI  section  arc  restored  to  the  program.  That  is.  Ci.  the  platform  torquer  gain,  is  5.0  It- 


14 


■MM 


lbs;  K.  the  inverse  of  the  platform  inertia  is  30.  (Ib-see  -ft  I I I R I and  I IR2  are  made  2 ft-lbs; 
the  spring  constant  CM I’l  and  the  viscous  damping  If  remain  at  12.  ft-lb/rad  and  0.1  It- 
1b  rad/sec,  respectively  I he  spring-mass-dumper  system  no  longer  resembles  a seeker 
platform  1 he  second  order  system  has  a damping  coefficient  of  about  0.08,  and  an  undamped 
natural  frequency  of  about  19  rad  sec.  I his  has  added  to  it  a healthy  but  not  overwhelming 
frictional  torque  of  2 ft-lbs.  IheRMN  parameter  for  model  B is  10  '"rad/ sec.  The  simulation 
is  run  by  invoking  GO.  and  the  resulting  angles  and  angular  rates  are  shown  in  Figures  21-23. 
Figure  24  compares  the  angles  for  models  A.  B and  C.  Again  one  notices  models  A and  B do 
not  reach  a steady  state  value  for  a step  input.  However,  their  values  are  close  due  to  MAXI 
being  set  at  2 5 msec  Model  It  w as  looked  at  with  R M N set  at  0.0025  and  0.05.  and  the  plots  are 
shown  in  Figures  23  and  2f>.  respectively.  The  importance  of  picking  RMN  for  good 
performance  of  model  B is  readily  apparent.  The  response  of  the  system  without  friction  is 
shown  in  Figure  27. 

K is  set  to  300  to  increase  the  undamped  natural  frequency  to  60  rad  sec.  B is  ntadeO.032  to 
keep  the  damping  coefficient  0.08.  Fven  though  the  integration  step  si/e  of  2.5  msec  is  still 
adequate  for  stability.  Figures  28-31  show  that  the  performance  of  Models  A and  B has 
degraded.  Figure  32  shows  the  response  of  the  system  with  no  friction. 

With  K = 30  and  B = 0.l  a 1.0  H/ sine  wave  is  used  as  the  system  forcing  function.  G remains 
5 ft-lb,  RMN  is  10  rad  / sec  and  MAX  T is  2.5  msec.  Figures  33-33  are  the  results  for  models 
C.  Band  A.  respectively.  Figure  36  shows  the  good  match  for  all  three  models.  Figure  3 7 shows 
the  response  for  the  system  without  friction. 

When  K is  made  300  and  B is  0.032,  the  accuracy  of  models  A and  B deteriorates  as  shown  in 
Figures  38-41.  Model  B had  RMN  set  at  10  Resetting  it  to  0.1  rad/sec  produced  the  nice 
match  with  model  C as  shown  in  Figures  42  and  43.  Figure  44  shows  the  system  response  for 
the  no  friction  case. 

With  k = 30..  B is  set  to  2.0  to  obtain  an  overdamped  system.  This  gives  a damping 
coefficient  of  1.58,  and  the  undamped  natural  frequency  remains  19  rad  sec.  The  frictional 
torque  remains  at  2 ft-lb.  Figures  43-47  show  the  results  for  a step  input  to  the  system  with 
models  t'.  B and  A.  respectively.  Figure  48  shows  the  good  overlays  obtained  for  all  three 
models.  This  is  due  to  the  angular  rate  never  dropping  below  zero.  Figure  4V  is  the  system 
response  without  friction. 

Setting  K to  300.  and  B to  0.632  makes  the  linear  undamped  natural  frequency  60  rad  sec. 
and  the  damping  ratio  remains  at  1 .58.  Figures  30-33  show  the  good  comparisons  for  the  three 
models.  Figure  34  shows  the  system  response  without  friction. 
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Although  models  A and  B have  angular  rates  with  high  oscillations  their  angle  histories 
compare  very  well  to  that  ol  model  ( as  show  n in  / inures  SS-5/i.  Figure  5V  is  the  frictionless 
response  of  the  system 

With  K = 300.  and  B = 0.632.  one  secs  significant  errors  in  models  A and  B as  shown  in 
/ igurex  M-fi.i  R MN  is  changed  tot).  I . and  the  resulting  improvement  in  model  B is  reflected  in 
Figures  <W  and  63.  Figure  lift  shows  the  system’s  frictionless  response. 

STOP 

terminates  the  simulation.  The  PRINT  file,  now  attached  to  the  interactive  terminal  as  a local 
lile,  is  hatched  for  line  printing  to  interactive  terminal  30. 

BATCH.  PRINT.  PRINT.  30.  HDMDD 

I he  Tektronix  terminal  is  logged  out  and  terminal  30  is  logged  in  to  obtain  the  line  printer 
listing  of  PRINT. 

Model  C is  clearly  superior  in  accuracy  to  models  A and  B under  the  conditions  tested 
However.  Model  B ,s  much  less  complex  and  gives  good  results  when  RMN  is  adjusted 
appropriately.  Angle  errors  in  A and  B with  a sine  wave  forcing  function  were  relatively 
smaller  than  those  tn  which  the  input  was  a step  except  for  the  overdamped  case.  Obviously  the 

greater  the  relative  magnitude  of  the  frictional  torque,  the  greater  the  accuracy  required  in  the 
modelling  ol  friction. 


CMPL 


Figure  1.  The  test  problem  used  to  study  the  friction  models. 


Figure  3.  Friction  Model  B implemented  with  two  submodels  that  are  switched  in 
and  out  as  a function  of  >!/. 


Figure  4.  Friction  Model  C. 
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Figures.  Listing  of  the  ACSL  FRICTN  macro  code  for  Model  C. 


Figure  6.  Flow  diagram  for  the  FRICTN  macro. 


PR  OV ICES  ENVIRONMENT  FOB  FRICTN  MODULE 


PROGRAM  FRICTION  TEST 


INITIAL 

LOGICAL 

CONSTANT 


CUMP 

OUMP=. FALSE.  • 

. TWOPI.fc.2A31  AS 
C INTERVAL  CINT  * 0.410 
ALGORITHM  IALG  « A 
NSTEPS  NSTP  * I 
MAXTEOVAL  MAX'?  . 0 . 402S 
LOGICAL  SINSTP 

CONSTANT  B « 0.1  . 

. TFR?  > 2.0  . 

. SIIC  ■ 0.4  . 

. G * S.O  . 

. SIPPIC  * 4.0  . 
WP  « M* TMOP I 


END 


Rmn.I.E-34 


CmPl  * 1?.0 
a « 30.0 
ICP  * 0.0 
m * 1,0 
ICPP  « 4.0 


. ISTP  . 2.0 


. TFRi  * 2.0 
. IC  = 0.0 
, SIPIC  * 0.0 
, SINSTP  . .FALSE. 


DYNAMIC 

DERIVATIVE 

.. STEP  COMMANO  TORQUE  IN 

TM3  « G«RSW(SINSTP.  SIn<*P*T).  STEP(O.l)) 

" SUM  FOR  NET  T0B0U2  MINUS  FRICTION 

TNET  ■ TM3  - B«SID  - CMpL«SI 

m OUTPUT  RATE  ♦ ANGLE  FROM  FRICTN  MACRO 

FRICTNIS10.  TNET.  K.  TER  1 » TFR2.  IC> 

SI  « INTEGiSID.  SIIC) 


output  rate  ♦ angle  from  another  model 

TNFT3  ■ TM3  - M*SIPC  - CMPL»SIP 

TNETP  ■ RSWIAPSISIPO.LT.RMN'  OEADI-TFRI.  TFRI,  TNET2)  * TNET2 
♦ SIGNITER2.-SIP0) ) 

SIPC  ■ INTfG  < TNETP •«  « ICPl 

SIP  « INTJGISIPO.  SIPIC) 


OUTPUT  RATE  ♦ ANGLE  FROM  AN  OFTEN  USEO 

FRICTION  MODEL 

TNET3  ■ TM3  - 8«SIPP0  - CMPL«SIPP  ♦ SIGN (TFR?. -5IPPD) 

SIPPO  ■ INTEG (K»TMET  3.  ICPP) 

S1PP  « INTEG (SIPPO.  SIPPIC) 


END 

TERMTIT  ,GT . TSTP) 

ENO 

TERMINAL 

IF(OUMP)  CALL  OEBUG 

END 

ENO 

0000004000004044004004 

0000000000000044004000 


Figure  7.  Listing  of  the  ACSL  test  program. 


CONN*-' C T . Ai'TwuT 
FTL.WC 

ATTACH*  INP'iT  .TF*  • ir- 

ATTACH. i <~,au  , TfK,  in  = n(')f  < »n  . y = 1 

l r, 0° 


SC  T TF0  1sl4..TF«3=U,.T;r->n.  ,<  = 1 0 . . MMT  =0 . 0 1 h 
SFT  T ! Tl  F ( ">)  ="M  vM  IN|  7h" 
of  o T "*H|"a?,,si.Mn 
Pt  OT  MP.MP  1 
PLOT  SlPP.MPot' 

PLOT  STPO.STP.SI 


STOP 

HATCH.PPIM  .ppI'T  , Hf  *jr  • H[)Mnp 


Figure  8.  Outline  of  the  Tektronix  keyboard  Inputs  tor  interactive 
session  setup  and  ACSL  run  time  commands. 


oooooooooooooooooooooe 

HONDO. 077000, 

ACCT  • • • 

ATTACH <MACFIL.DCMACFlL*ID*DCACSLSYS> 
ATTACH (ACSL .DCACSL. I0=0CACSLSYS) 

ACSL  ( IMNPUT) 

RETURN. ACSL. MACFIL. 

FTN(I=C0MPILE.R=2> 

MAP. OFF. 

REOUEST.LGOB. «PF. 

ATTACH (ACSLL IB. OCACSLL IB *I0=DCACSLSYS> 
ATTACH,PLT.TCKTRONIXAOH.ID*WTPLOT.CY=2. 
LOSET(lIB*PLT.SUBST=ZZORAW-TEKPLT» 
LOSET(LIB=ACSLLIB.PRESET=INOEF) 

LOAD.LGO. 

NOGO.LGOR. 

RETURN .ACSLL IB .PLT • 
CATAL0G.LG0B»TEK.ID=DCXXXM*CY*1. 

EXIT. 

0000000000000000000000 


Figure  9a.  Job  control  stream  used  to  generate  LGOB,  the  absolute 
binary  file  of  the  test  program. 


SET  PRN«9 

SET  prnplt=.f4*calplt*.t.,ttlcpl=.t. 

SET  TITLE»"FRICTN  TESTER 
SET  OUMP».T. 

SAVE  "BASE" 

PROCEO  GO 
START 

PRINT  "ALL" 

ENO 

PREPAR  T.SI.SI0.SIP.SIPD.SIPP*SIPP0.TH3,TNET,TNET?.TNET3 
SET  CM0*0IS 


Figure  9b.  ACSL  run  time  commands  on  file  INPUT. 
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Figure  10.  Angle  and  rate  response  ot  a simulated  platform  gimbal  to  a step  input  — model  C. 
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response  of  a simulated  platform  gimbal  to  a step  input  — model  A. 
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Comparison  ol  angle  responses  of  a simulated  platform  gini' 
bal  to  a step  Input  for  models  A,  B and  C. 
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response  of  a simulated  platform  glmbal  to 
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Figure  14b.  Model  B and  C angle  response  comparison  (or  simulated 
platform  glmbal  to  a step  input  — RMN  = 0.015. 
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response  of  a simulated  platform  gimbal  to  a 1.0  Hz 


FRICTH  tester  - 
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response  of  an  underdamped  second  order  system  to  a step  input  — model  C. 
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response  of  an  underdamped  second  order 
tp  input  — model  B. 
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Comparison  of  angle  responses  for  an  underdamped  second 
order  system  to  a step  input  for  models  A,  B and  C. 


RUN  3.05000000 
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Figure  27.  Angle  and  rate  response  of  the  frictionless  underdamped 
ond  order  system  to  a step  input. 
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response  of  a high  gain  underdamped 
o a step  input  — model  B. 
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response  of  a high  gain  underdamped 
> a step  input  — model  A. 
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Comparison  of  angle  responses  of  a high  gain  underdamped 
second  order  system  to  a step  input  for  models  A,  B and  C. 
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Figure  32.  Angle  and  rate  response  of  the  frictionless  high  gain  under 
damped  second  order  system  to  a step  Input. 
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•e  33.  Angle  and  rate  response  of  an  underdamped  second  order 
system  to  a 1 .0  Hz  sine  wave  — model  C. 
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response  of  an  underdamped  second  order 
Hz  sine  wave  — model  B. 
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Comparison  of  angle  responses  of  an  underdamped  secon 
order  system  to  a 1 .0  Hz  sine  wave  for  models  A,  B and  C. 
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Figure  37.  Angle  and  rate  response  of  a frictionless  underdamped 
ond  order  system  to  a 1.0  Hz  sine  wave. 
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Figure  38.  Angle  and  rate  response  of  a high  gain  underdamped  second 
order  system  to  a 1.0  Hz  sine  wave  — model  C. 
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response  of  a high  gain  underdamped  second 
1.0  Hz  sine  wave  — model  B. 
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Figure  40.  Angle  and  rate  response  of  a high  gain  underdamped  second 
order  system  to  a 1.0  Hz  sine  wave  — model  A. 
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Comparison  of  angle  responses  of  a high  gain  underdamped 
second  order  system  to  a 1.0  Hz  sine  wave  for  models  A,  B and  C. 
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Figure  42.  Angle  and  rate  response  of  a high  gain  underdamped  second 
order  system  to  a 1.0  Hz  sine  wave  for  model  B with  RMN  = 0 
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Figure  44.  Angle  and  rate  response  of  a frictionless  high  gain  under- 
damped second  order  system  to  a 1.0  Hz  sine  wave. 
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Figure  45.  Angie  and  rate  response  of  an  overdamped  second  order  sys- 
tem to  a step  input  — model  C. 
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response  of  an  overdamped  second  order  sys' 
t — model  B. 
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Figure  47.  Angle  and  rale  response  of  an  overdamped  second  order  sys 
tem  to  a step  input  — model  A. 
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B and  C. 
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Figure  50.  Angle  and  rate  response  of  a high  gain  overdamped  second 
order  system  to  a step  input  — model  C. 
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response  of  a high  gain  overdamped  second 
tep  input  — model  B. 
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Figure  52.  Angle  and  rate  response  of  a high  gain  overdamped  second 
order  system  to  a step  input  — model  A. 
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responses  of  a high  gain  overdamped 
to  a step  input  for  models  A,  B and  C. 
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igure  54.  Angle  and  rate  response  of  a high  gain  friction!^ overdamped  second 
order  system  to  a step  input. 
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response  of  an  overdamped  second  order  sys 
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response  of  an  overdamped  second  order  sys 
«e  wave  — model  A. 
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order  system  to  a 1.0  Hz  sine  wave  — model  C. 
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Figure  61.  Angle  and  rate  response  of  a high  gain  overdamped  second 
order  system  to  a 1.0  Hz  sine  wave  — model  B. 


ERICTN  FESTER  - 8 JUNE  79 


Figure  62.  Angle  and  rate  response  of  a high  gain  overdamped  second 
order  system  to  a 1.0  Hz  sine  wave  — model  A. 
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responses  of  a high  gain  overdamped 
i to  a 1.0  Hz  sine  wave  for  models  A,  B and  C. 
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Figure  66.  Angle  and  rate  response  of  a frictionless  high  gain  over- 
damped second  order  system  to  a 1.0  Hz  sine  wave. 
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